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Abstract. Fundamental diagrams of vehicular traffic flow are generally multi- valued in 
the congested flow regime. We show that such set-valued fundamental diagrams can be 
constructed systematically from simple second order macroscopic traffic models, such as 
the classical Payne- Whitham model or the inhomogeneous Aw-Rascle-Zhang model. These 
second order models possess nonlinear traveling wave solutions, called jamitons, and the 
multi-valued parts in the fundamental diagram correspond precisely to jamiton-dominated 
solutions. This study shows that transitions from function-valued to set-valued parts in 
a fundamental diagram arise naturally as intrinsic properties of well-known second order 
models. 



1. Introduction 

In this paper a connection between set-valued fundamental diagrams of traffic flow and 
traveling wave solutions in second-order inviscid traffic models is presented. 

1.1. Types of Traffic Models. A wide variety of mathematical models for vehicular traffic 
exist. Microscopic models (e.g., [45t I39j ) describe the individual vehicles and their interac- 
tions by ordinary differential equations. Mesoscopic models (e.g., ^231 HU |29l E]) employ 
a statistical mechanics perspective in which vehicle interactions are modeled in analogy to 
kinetic theory Macroscopic models (e.g., [3aiMllinillll[MlE31E3[Ml[2|,as well as the dis- 
cussion below) describe traffic flow using suitable adaptations of the methods of continuum 
mechanics [Mj, which yields equations mathematically similar to those from fluid dynamics. 
There also exist other types of models that are not necessarily based on differential equations, 
such as probabilistic models (e.g., [Ij), and cellular models (e.g., [37]), which divide the road 
into segments and prescribe laws for the propagation of vehicles between cells. A detailed 
overview of the various types of models is given in f22J. 

All of these types of traffic models are related. For instance, mesoscopic models can be 
derived from cellular models [37] with stochastic components [1]; macroscopic models can 
then be interpreted as limits of mesoscopic models |38l [29] ; cellular models can be interpreted 
as discretizations ("cell transmission models") of macroscopic models [8, ,9j and thus relate to 
kinematic waves |llj ; and microscopic models are equivalent to discretizations of macroscopic 
models in Lagrangian variables. Through these relations, the macroscopic modeling results 
presented in this work may have consequences for various other types of traffic models. 

In this paper, we focus solely on inviscid (aggregated) single-lane macroscopic models for 
traffic flow on unidirectional roads. These types of models were started with the work of 
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Lighthill & Whitham [35J and Richards [36]. The Lighthill-Whitham- Richards (LWR) model 
describes the temporal evolution of the vehicular density along the road by a scalar hyperbolic 
conservation law, and is therefore denoted a "first order model". The first "second order 
model", due to Payne [l2] and Whitham [5l], introduces the vehicular velocity field as an 
independent variable. As a reaction to certain worrisome mathematical properties of the PW 
model [lOj, Aw & Rascle [2J, and independently Zhang [55j devised a different second order 
model, which has been subsequently extended and generalized [T8t [3]. 

In all of these models, traffic is described by inviscid conservation or balance laws. There 
also exist macroscopic models with higher spatial derivatives: KdV-type models (e.g., |3 11130] ) 
involve a dispersive term, and viscous models (e.g., [271 [28]) possess a viscosity term. These 
two classes of models explicitly describe the dynamics of the traffic flow on multiple length 
scales and avoid discontinuous solutions. In contrast, inviscid models replace thin transition 
zones (e.g., upstream ends of traffic jams) by shocks, and the specific dynamics inside the 
transition zone are replaced by appropriate Rankine-Hugoniot jump conditions [12j . 

1.2. Fundamental Diagram. The fundamental diagram (FD) of traffic flow is the plot of 
the flow rate (i.e., the number of vehicles passing a specific position on the road per unit 
time) versus the vehicle density (i.e., the average number of vehicles per unit length). The 
first FD measurements were done in the 1930s by Greenshields [20]; they led him to propose 
a quadratic function Q{p) = pumax(l — p//Omax) describing the relation between fiow rate 
and density. Here, Umax is the velocity that a driver would assume when alone on the road, 
and /Omax is a maximum density at which the traffic fiow stagnates (bumper-to-bumper plus 
a safety distance). This form of the FD was later also used by Richards [IS]. However, 
subsequent measurements, such as those of Wardrop & Charlesworth |52j and H. Greenberg 
|17j . illuminated that a quadratic form is not an accurate description of true FD data, since, 
for instance, the maximum fiow rate generally occurs at densities significantly below ^Pmax- 
In reaction, different types of FD functions were proposed, such as by Lighthill &: Whitham 
[55] . Greenberg [lYj, Underwood [49j, Newell ^0], Daganzo |8], and Wang & Papageorgiou 

However, a more elementary shortcoming of function-based FD models is that true fiow 
rate vs. density measurements tend to exhibit a clear systematic spread: for sufficiently large 
densities, one value of p corresponds to multiple values of Q. The described "spreading" effect 
can be seen in Fig. [T| which displays a plot of flow rate vs. density, obtained from single loop 
sensor measurements at a flxed position on a highway. The data is taken from the RTMC 
data set [41] of 2003, provided by the Minnesota Department of Transportation (Mn/DOT), 
from a sensor on the southbound direction of I-35W (Minneapolis, MN). Each point in the 
diagram corresponds to a time interval of At = 30s. The recorded flow rate equals the number 
of vehicles that have passed over the sensor during the time interval (the selection of discrete 
flow rate values is due to the integer nature of the number of vehicles); and the density is 
obtained from the fraction of time that a vehicle has blocked the sensor. For the plot in Fig. [T| 
we assume that the average vehicle length is 5m, and that the possible maximum density is 
achieved at 1/7. 5m, i.e., a safety distance of half a vehicle length is the absolute minimum 
that drivers will accept. 

Motivated by the spread of data in the FD, traffic fiow is often classified into regimes or 
phases. The two-phase traffic theory (see, e.g., [22l |50j) divided traffic flow into free flow for 
low densities, and congested flow for large densities. Three-phase traffic theory, originally 
developed by Kerner [25], distinguishes free flow for low densities, synchronized flow for 
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Figure 1. Fundamental diagram, obtained from 
sensor measurement data, aggregated over time 
intervals of At = 30s. The horizontal structures 
are due to the integer nature of the flow rate data. 
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Figure 2. Fundamental phase diagram, as fre- 
quently depicted in the traffic literature (e.g., 
PSI [5U]). Shown here is the free flow curve, 
and a synchronized flow region. 



medium densities, and wide moving jams for high densities. A possible depiction of traffic 
phases for the sensor data is shown in Fig. [2| The free flow region can be characterized by 
a relatively well defined single valued curve (with the spread due to data noise), while for 
medium densities (synchronized flow) a multi- valued flow region appears to be needed. 

By design, first order traffic models assume a function-based FD, i.e., each density value 
p corresponds to a unique fiow rate value Q. Therefore, more general traffic models are 
required to explain the observed spread in FDs. In recent years, a variety of approaches has 
been presented that describe the spread in FDs by set-valued traffic models. One class of 
approaches are phase transition models [Si [3 Sj that explicitly prescribe one traffic model 
for low densities, and another traffic model for larger densities, and construct mathematical 
transition conditions that guarantee well-posed solutions. Other models are based on inserting 
an instability into a second-order model at densities that are known to lead to instability in 
reality ^47]. 

1.3. Contribution of this Paper. In this paper, we demonstrate that the observed struc- 
ture of FDs, namely a unique Q vs. p relationship for small densities and a set-valued relation- 
ship for larger densities (see Fig. [T]), does already arise within classical second order models, 
and thus does not necessarily require explicitly prescribed phase transitions. Specifically, we 
shall develop the theory for the classical Payne- Whitham model |42l I54j . and then extend 
the results to the more recent inhomogeneous Aw-Rascle-Zhang model [2| 1551 [T8] . The key 
observation is that second order models with a relaxation term possess an intrinsic phase- 
transition-like behavior: for certain densities, the solutions are very close to the solutions of a 
first order model (as made precise in the introduction to Sect.js]), while for other densities, the 
solutions develop jamitons [15], i.e., self-sustained, attracting, traveling wave solutions with 
an embedded shock (described in more detail in Sect.|4]). As we shall show, these traffic waves 
result from instabilities and thus can be interpreted as manifestations phantom traffic jams; 
and the resulting jamitons can be regarded as models for stop-and-go-waves. Both of these 
phenomena are observed in real traffic flow [21], and have been reproduced experimentally 
jUj. That being said, we emphasize that the key statement of this paper is that set-valued 
FDs are already built into second order models of traffic flow with a relaxation term, such 
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as the PW model and the inhomogeneous ARZ model. No statement is made about whether 
these models are better or worse than other (possibly more complicated) models in describing 
real vehicular traffic flow. 

This paper is organized as follows. Sect. [2] consists of a brief presentation of the two types 
of second order traffic models studied here, as well as their mathematical properties that are 
relevant for the construction of jamitons. In Sect. [Sj we outline a key theorem in the study 
of second order systems of hyperbolic balance laws with a relaxation term: the equivalence 
of the instability of uniform base states, and the sub-characteristic condition. Then, for 
common choices of parameters in the considered second order traffic models, the implications 
of this relationship are outlined. The explicit construction of jamiton solutions in the two 
types of second order traffic models is presented in Sect. |4j In particular, the equivalence 
between the breakdown of the sub-characteristic condition and the existence of jamitons is 
shown. In Sect.jSj we demonstrate how the preceding results about jamitons and the stability 
of uniform traffic flow can be employed to systematically construct set-valued fundamental 
diagrams that qualitatively reproduce the key features of fundamental diagrams obtained from 
measurements. Finally, in Sect. [6j conclusions are drawn and an outlook towards extensions 
of the presented work is given. 



2. Second Order Traffic Models with Relaxation Terms 

Macroscopic traffic models describe the temporal evolution of the vehicle density p{x,t), 
where x is the position along a road (multiple lanes are aggregated into a single variable) , and 
t is time. Letting u{x, t) denote the mean flow velocity of the moving vehicles, the density 
function evolves according to the continuity equation 

Pt + {pu)x = . (1) 

First order traffic models, such as the Lighthill-Whitham-Richards (LWR) model [35l H6] . 
assume a functional relationship between p and u, i.e. u = U{p). We call U the desired 
velocity function, and assume it to be decreasing with p. A simple choice is 

U{p) = nmax(l - P/Pmax) , (2) 

which is a linear interpolation between a maximum velocity limax of a single vehicle on an 
otherwise empty road (i.e., p = 0) and a stagnated flow at some maximum density, i.e., 
U{pniax) = 0. Because Q does not match very well with measurement data (see Fig. [T]), a 
variety of other desired velocity functions has been proposed [T71 HHl 1301 [H [51] . The choice 
of a desired velocity function U{p) leads to a scalar hyperbolic conservation law model 

Pt + {Q{p)). = , (3) 

where the flux function Q{p) = pU{p) describes a unique Q vs. p relation in the associated 
fundamental diagram of traffic flow. In contrast, second order models augment ([T]) by an 
evolution equation for the velocity held u. Below, we outline two classical types of second 
order models: the Payne- Whitham model, and the inhomogeneous Aw-Rascle-Zhang model. 
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2.1. Payne- Whitham Model. The Payne- Whitham (PW) model [121 [5l], written in the 
primitive variables p and n, reads as 

pt + {pu)x = , 

(4) 

ut + nn^ + -pP{p)x = 7 {U{p) - u) , 

where the traffic pressure p{p) is a strictly increasing function and r is a relaxation time that 
dictates how rapidly drivers adjust to their desired velocity U{p). The PW model Q is a 
system of hyperbolic balance laws that possesses a relaxation term in the velocity equation. 
In analogy to fluid dynamics, the PW model is generally taken to possess the conservative 
quantities p and q = pu, in which it reads as 

Pt + qx = , 

(5) 



qt + {'^+p{p))^ = HpU{p)-q) 



While Q and ([s]) are equivalent for smooth solutions, the propagation of shocks is determined 
by the Rankine-Hugoniot conditions [12j associated with the form ([s]). The characteristic 
velocities of the PW model are 

Ai = 2-v9(p) and A2 = ^ + v9(p), (6) 

i.e., in the PW model, information travels with the velocities u it c{p), where c{p) = \/p'{p)- 

System ([s]) is the PW model written in Eulerian variables p{x, t) and q{x, t). An alternative 
description, as for instance employed in [19], is to use Lagrangian variables v{a, t) and u[a, t). 
Here a is the vehicle number, defined so that 

da = pdx — pudt , (7) 

and V = 1/p is the specific traffic volume, i.e., the road length per vehicle. In these variables, 
the PW model reads as 

Vf - = , 

(8) 



ut + {p{y)). = ^{u{y)-u) 



where j5(v) = p(l/v) and U{v) = U{l/v). For notational economy, we shall hereafter omit 
the hat in these two functions, except when required for clarity. 

The Rankine-Hugoniot jump conditions associated with ([s]) are 

m [vl — lu] = , 

(9) 

m [u\ + [p{'v)\ = . 

Here [z] is the magnitude of the jump in the variable z across a shock, and —m is the 
propagation speed of the shock (in the Lagrangian, i.e., vehicle-centered, frame of reference). 

2.2. Inhomogeneous Aw-Rascle-Zhang Model. The inhomogeneous Aw- Rascle- Zhang 
model ^ i55t il8j in primitive Eulerian variables reads as 

Pt + {pu)x = , 

(n + h{p))t + u{u + h{p))x = \ {U{p) - u) , ^^^^ 
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where the hesitation function h{p) is a strictly increasing function. As in the PW model, 
one assumes, in a loose analogy to fluid dynamics, the conservative quantities to be p and 
q = p{u + h{p)). In these variables, the model reads as 

Pt + qx = , 

Qt+{i- h{p)q)^ = \ (piUip) + h{p)) - q) . ^^^^ 

The characteristic velocities of the ARZ model are 

Ai = 2 - h{p) - ph'{p) = u- ph'{p) and A2 = ^ - h{p) = u , (12) 

where u = ^ — h{p). Thus, in the ARZ model, information propagates along one family of 
characteristics with the vehicles, and along the second family of characteristics against traffic 
(relative to the vehicles) with velocity ph'{p). 

In Lagrangian variables, the ARZ model becomes 

(13) 



(v + Mv))i = i(?7(v)-n) 



where h{v) = h{l/v) and U{v) = U{l/v). Again, we shall generally omit these hats from now 
on. The associated Rankine-Hugoniot conditions are 

m [vl — lu] = , 

(14) 

[u] + [h{v)] = . ^ > 

2.3. Assumptions on the Desired Velocity, Pressure, and Hesitation. In addition 
to requiring that U{p), p{p), and h{p) be twice continuously differentiable, we impose the 
following assumptions: 

(a) > and < 0. This is equivalent to U{p) being decreasing, and the flux function 
Q{p) = pU (p) in the corresponding LWR model Q being concave. The concavity of 
Q{p) is a frequent assumption in traffic flow (see for instance the discussion in [35]). In 
particular, it guarantees a nice shock theory for the LWR model, where the shocks satisfy 
the Lax entropy condition [32] strictly. 

(b) -gf < and > 0. These conditions lead to well-posed equations for the PW model 
m^, even in the presence of shocks [32]. These conditions are equivalent to: as functions 
of p, the pressure p{p) is increasing and the function pp{p) is convex. 

(c) -g^ < and > 0, for the same reasons as for p. 

Note that all the example functions U, p, and h used in this paper satisfy these assumptions. 



3. Instability of Uniform Flow and Sub-Characteristic Condition 

A particular feature of balance laws of the form Q and (10) is the relaxation term in the 
(generalized) momentum equation, which is expected to drive the actual velocity u towards 
the desired velocity U{p) when r ^ 1. If this is indeed the case, the continuity equation ([l]) 
shows that the solution will be very close to the one of the LWR model ([s]) with the flux 
function Q{p) = pU{p). 
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Note that ^ has a single characteristic velocity 

^^ = Q'{p) = Uip) + pU'{p), (15) 
while the second order models Q or (10) possess two characteristic velocities Ai < A2, which 



depend in particular on the traffic pressure p{p), or the hesitation function h{p), respectively. 
Moreover, one can observe that all three models Q, (10), and ([s]) share a set of solutions, 



namely the base state solutions, in which p{x,t) = p = const, and u = U{p) = const. These 
solutions represent uniform traffic flow with equi-spaced vehicles (whose midpoints are a 
distance 1/p apart). Let now p = p, u = u = U{p) be any given base state solution of a 
second order model and consider the following two conditions that can apply: 

LS (linear stability) 

The base state (/>, u) is linearly stable, i.e., any infinitesimal perturbations to it decay 
in time. 

sec (sub-characteristic condition) 

The characteristic velocity p = p{p), given by ([15|, of the reduced LWR model ^ 



falls in between the two characteristic velocities Xj = Xj{p,u) of the second order 
model, i.e. 

Ai < /i < A2 . (16) 
These two conditions are related by the following 

Theorem 1 (Whitham). The two conditions LS and SCC are equivalent. 



The proof, given by Whitham in 1959 [53], and thoroughly reviewed in his 1974 book 
Chap. 10], is constructive: for a general 2x2 hyperbolic conservation law with a relaxation 
term in one equation, the standard linear stability analysis (calculating the growth factor 
for basic wave perturbations of the form e*''^) is performed, and subsequently related to the 
sub-characteristic condition. 

Besides the theoretical insight that Thm. [l] yields (see the discussion in [54j ) , it allows one 
to substitute the task of performing a linear stability analysis by, for instance, normal modes 
(e.g., as done for (j4]) in [T^l App. B] and for ( [To| ) in [TO]) by checking the sub-characteristic 
condition, which is often times simpler and faster to verify. 

In the context of traffic flow, an important question is: when can one use a first order model 
and (essentially) get the same answer as that given by a second order model? This question, 
which turns out to be intimately related to the result in Thm. [T| has been investigated by 
several authors, e.g., Liu [36]) Chen, Levermore & Liu [5], and Li & Liu [3l] (note that these 
papers consider situations that are slightly more general than the ones arising for traffic flow) . 
Thus, for example, it is known that: 

1. Solutions of a second order traffic model with uniform base state initial conditions, plus 
a sufficiently small perturbation, converge towards this base state (i.e., the perturbation 
decays) as t — )■ 00, if the base state satisfies the SCC. In other words, the equivalence in 
Thm. [Tjcan be extended from relating SCC and LS, to relating SCC and stability. 

2. The solutions to second order models that satisfy the SCC everywhere, converge towards 
solutions of the reduced LWR model in the limit r — )■ 0. More generally, if r is not 
small, then the solutions converge towards the solutions of a reduced LWR model with a 
nonlinear "viscosity", which is of magnitude 0{t). 

An important question that the work cited above does not address is: 
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3. If constant base state (plus small perturbation) initial conditions (strictly) violate the 
sec, how do the solutions to a second-order model evolve in time? 

To the authors' knowledge, no rigorous results are known in this case, except for the fact that 
the solutions to the second order models need not look like solutions to the reduced LWR 
model even if r ^ 1. However: 

4. Numerical studies (for the PW model in [15], and for the ARZ model in [19j) indicate 
that the solutions of second order traffic flow models tend towards a jamiton-dominated 
regime, characterized by the presence of self-sustained finite amplitude traveling waves. 

The construction of these traveling waves is addressed in Sect. |4j However, before doing so, 
we calculate the sub-characteristic condition for the PW and the ARZ model below. 

Remark 1. The situation for the degenerate case of the SCC condition, i.e., Ai = /U < A2 or 



Ai < /i = A2, in (16), seems to be considerably more involved and we do not consider it in 
this paper. 

3.1. Payne- Whitham Model. Using the characteristic velocities of the PW model, 
and of the LWR model, (15), the sub-characteristic condition (16) reduces to 



-vV(p) < pu'ip) < ^Mp) . 

Because U{p) is decreasing and p{p) is increasing, the second inequality is always satisfied, 
rendering the SCC, and thus the stability of uniform flow of density p, equivalent to 

> {U'{p)f . (17) 

Regarding reasonable choices for U{p) and p{p), the desired velocity function has a clear 
physical interpretation. In this paper, we consider two forms of U{p): the simple linear 
relationship ([2]), and a more complicated shape, given via (33). In contrast, which choices 



are reasonable for the traffic pressure is more difficult to motivate. For the linear velocity 
U[p) = Umax(l — p//Omax), we briefly study two classes of traffic pressures in terms of their 
behavior with respect to the (in)stability of base states. We argue that reasonable traffic 
models for phantom traffic jams must have the property that low density base states are 
stable {free flow), instability occurs above a certain critical threshold density [congested flow 
or synchronized flow), and possibly stability re-occurs for densities close to the maximum 
density pmax {wide moving jams) , where the terminology is in line with [221 EH]. We distinguish 
two types of pressures: 

(A) regular pressure: Drawing an analogy to gas dynamics, a simple choice is p{p) = f3p'^ , 

— "^ (u )^ 

where 7 > 0. In this case, the SCC becomes p^-^ > -^T^- One can see that for 
7 > 3, low densities violate the SCC, and for 7 = 3, no transitions between stability or 
instability occur. Hence, reasonable behavior occurs only for 7 < 3, for which instability 

occurs above the critical density pcvit = (/37(Pmax)^/('Umax)^) • Some possible choices 
are 7 = 1, as in [26], or 7 = 2, which renders the PW model ^ equivalent to the shallow 
water equations with a relaxation term. 

(B) singular pressure: Stability for densities close to pmax can be achieved by selecting a 
pressure p{p) with a pole at pma.x- For example, as in [T5j, we can take 

P{P) = -/5(P/Pmax + log(l - p/Pmax)) , (18) 
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for which the SCC becomes , > or equivalently ^(1 ^) < ^ — ^ 

P(Pmax-pj (Pmax) ' ^ Pmax ^ Pmax ' Pmax(Mmax) 

Thus, there is a range of unstable base states that is centered around ^Pmux, and is of 
width ((ipmax)^ " ^Pmax \ ^ ^ Hence, for low densities, and high densities (i.e, close to 
Pmax) 5 uniform traffic flow is stable. Thus, singular pressures can generate traffic regimes 



qualitatively close to those presented by Kerner [25J. We shall use this pressure (18) for 
some of the examples presented in Sect. [5] 

3.2. Inhomogeneous Aw-Rascle-Zhang Model. Using the characteristic velocities of the 



ARZ model, (12), and of the LWR model, (15), the sub-characteristic condition (16) reduces 
to 

-ph'ip) < pU'ip) < . 

Because U{p) is decreasing, the second inequality is always satisfied, rendering the SCC, and 
thus the stability of uniform flow of density p, equivalent to 

h'ip) > -U'ip) . (19) 

As in the PW model, we study two types of hesitation functions, in conjunction with the 
desired velocity function U{p) = Wmax(l — p/Pmax): 

(A) regular hesitation: For the choice h{p) = f3p'^ with 7 > 0, the SCC becomes p'^~^ > 
■g^^^^. Similarly as with the PW model, realistic behavior of (in)stability regimes only 



occurs if 7 < 1. 

(B) singular hesitation: As with the PW model, stability can be recovered close to pmax 
by letting the h{p) become singular there. A possible choice, suggested in [3], is 

hiP) = ^{T^y- (20) 

Since for p ^ 1, the function h{p) behaves as in the regular case, one can again only 
achieve realistic (in)stability regimes if 7 < 1. In this case, h{p) possesses an inflection 
point and thus can generate the traffic regimes presented by Kerner [25j. In Sect. [5| we 
shall employ ( 20 ) with 7 = | • 



4. CONSTRUCTION OF JAMITON SOLUTIONS 

Jamitons are self-sustained traveling waves in second order traffic flow models with a re- 
laxation term. They possess the same mathematical structure as detonation waves in the 
Zel'dovich-von Neumann-Doring (ZND) theory [TJ] of reacting gases. A ZND detonation 
wave consists of a compressive shock with a smooth reaction zone attached to it, which is 
traveling at the same velocity as the shock. The ideas presented in this section are an adapta- 
tion of this theory to the PW and the ARZ traffic models. A theory valid for generic second 
order traffic flow models is presented in a companion paper |16j . 

In order to obtain expressions for jamitons, a traveling wave ansatz is performed. This 
approach has been presented in |15j in Eulerian variables, wherein one considers solutions of 
the p{x,t) = p{t]) and u{x,t) = u{r]), where rj = is the self-similar variable. Here, we 
use the Lagrangian variables defined in Sect. [2j which have in particular been employed in 
Greenberg's study of the ARZ model |I9] . 



10 



B. SEIBOLD, M. R. FLYNN, A. R. KASIMOV, AND R. R. ROSALES 



4.1. Payne- Whitham Model. We consider system ([s]) and seek solutions of the form 
v(x,t) = v(x) and u{x,t) = where x = ^2ii+£. jg the self-similar variable. For sim- 

plicity of notation, we shall omit the hats in the following. Substituting this self-similar form 
into ([8|) yields the relations 



'ru'ix) + p'(v(x))^'(x) = ^ (f/(v(x)) - n(x)) 



Relation (21) implies that 



mv 



u 



(21) 
(22) 

(23) 



where s is a constant of integration. 



Remark 2. In Lagrangian variables, —m is the propagation velocity of jamitons, i.e., jamitons 
travel backward at speed m with respect to the moving vehicles. The constant s has the 
interpretation of a "Lagrangian mass flux" . In Eulerian variables, the meaning of these two 
constants is exactly reversed: jamitons move at velocity s with respect to a fixed observer, 
and m is the "mass" flux of vehicles past a stationary observer. Note that due to the entropy 



condition described in Sect. 4.1.1, solutions of the type above can exist only for m > 0. On 



the other hand, there is no sign restriction on s. 



Using relation (23) to substitute u by v in (22 ), we obtain a first order ODE for v(x), which 
reads as 

^^(v(x)) 



v'(x) 



(24) 



r'(v(x)) ' 

where the two functions w and r are defined as 

w{'v) = U{y) — (mv + s) and r(v) = ^(v) + m^v , (25) 

and depicted in Fig. |4] and Fig. [3j respectively. At first glance, it appears that the constants 
m and s can be chosen independently. However, this is not the case. To see this, we use 



the assumptions p'(v) < and p"(y) > (see Sect. 2.3). Then, the denominator in (24), 
r'(v), may have at most one root, vs, such that p'(vs) = —m?. The ODE (24) can only be 



integrated through vs, if the numerator in (24), w(v), has a (simple) root at vs as well, i.e. 

mvg + s = C/(vs) . (26) 
In ZND detonation theory, relation (|26|) is known as the Chapman- J ouguet condition |14] . 



and the point at which r'(v) vanishes is called sonic point. We can therefore parameterize 



the smooth traveling wave solutions (of the jamiton ODE (24)) by the sonic value vg. Given 
vg, the two constants are 

= \/— p'(vs) and s = U (vg) — mvg . (27) 
Hence, we obtain a one-parameter family of smooth traveling wave solutions, parameterized 



by Vg, each of which is a solution of (24) 



Whether, in the vicinity of vg, the function v(x) is increasing or decreasing, is determined 
by the sign of the right hand side in ([24]), which, by L'Hopital's rule, equals 

^«'(vg) [/'(vg) 



m 



r"(vg) p"(vg) 

Because p"(v) > 0, we note that v(x) is increasing if C/'(vg) > m, and decreasing if C/'(vg) < 
m; the degenerate case C/'(vg) = m yields only constant solutions. 
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Figure 3. The r-v diagram in the p-v plane. 
The sonic value vs is located in the minimum of 
r(v), and the maximal jamiton (blue line) con- 
nects the states vr and vm- 
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Figure 4. The w-v diagram in the u-v plane. 

Jamitons exist exactly if w{v) possesses a second 
root Vm right of vg, which is the case if U'{vs) > 
m. 



Using that U"{v) < (see Sect. 2.3), we know about the numerator in (24), w{v), that: if 
U'{vs) > ni, a second root can exist for v > vs, but no further roots can exist for v < vg (see 
Fig. |4]). If, on the other hand, f/'(vg) < m, a second root can exist for v < vs, but no further 



roots can exist for v > vg. As a consequence, the solution of (24) can extend to infinity in 



the direction x — ^ +oo, but not into the direction x — ^ —oo (at a finite value of v). 

4.1.1. Shocks. Because smooth traveling waves cannot extend to infinity in both directions, 
they must contain a shock that connects a state (v~,u~) to the left of the jump to a state 
{v~^,u~^) to the right of the jump. A shock can be part of the traveling wave solution, if it 
travels with the same velocity m. The Rankine-Hugoniot jump conditions ^ now impose two 
restrictions on the states connected by the shock, as follows. The first condition in ^ implies 
that mv^ — u~ = mv^ — , i.e., the quantity mv — u is conserved across the shock. This is 



in line with condition ( 23 ) for the smooth part of a traveling wave. A second requirement is 
obtained by combining the two jump conditions in ^ to yield 

[v] + [p(v)] = , 

which is equivalent to the statement that r(v~) = r{v~^), where r is precisely the function 



defined in (25). Therefore, any jamiton shock connects two states v and v+ by a horizontal 
line in the r-v diagram, as shown in Fig. [5j 

In addition to the Rankine-Hugoniot conditions, a stable shock must also satisfy the Lax 
entropy conditions (see e.g. [MIE^)- For a shock that is part of a traveling wave (see [15j for 
a visualization), this requires that the specific volume decreases across a shock, i.e., v~ > v"*". 
Hence, in a smooth part of a traveling wave, v must be increasing with x^ ^iid therefore 
a necessary condition for the existence of jamitons is that C/'(vg) > m (see Fig. |4]). As a 
consequence, a jamiton solution v(x) corresponds to a smooth movement along the function 
r(v) in the p-v plane into the positive v direction, followed by a jump back from the right 
branch to the left branch along a horizontal line (see Fig. [5|. 

Remark 3. The PW model has been criticized [10| [2] for the fact that its solutions may 
involve shocks that travel faster than the vehicles. However, the arguments above point to 
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Figure 5. A jamiton in the r-v diagram in the 
p-Y plane. Below the maximal jamiton (blue line) 
lies an actual jamiton (brown line) that connects 
the states v+ and v^. 
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Figure 6. Jamiton function v(?7). An actual 
jamiton (brown curve) between v"'' and v~ lies 
on the maximal jamiton curve (blue), which con- 
nects vr and vm, the latter being at 77 — > oo. 



the following interesting fact: while such unrealistic shocks can certainly be triggered by 
particular choices of initial conditions, they cannot occur as parts of traveling wave solutions. 



4.1.2. Jamiton Construction. Based on the theoretical presentation above, next we present 
a systematic process for the construction of jamiton solutions, with an intuitive graphical 
interpretation. The construction is as follows: 

1. Select a sonic specific volume vs. 

2. Evaluate m = y/—p'{vs) and s = U{vs) — mvs- 

3. Draw the function w{'v) = U{y) — (mv + s) in the n-v plane (see Fig. |4|. 

4. If w'{vs) ^ 0, then no jamitons exist with this sonic value vs. Otherwise, find the unique 
root Vm > vs such that w{vm) = 0. 

5. Draw the function r(v) = p(v)-|-m^v in the p-v plane (see Fig.js]), and define rmin = r(vs) 
and Tmax = 'r(vM)- 

6. Find the unique value vr < vs that satisfies r(vR) = rmax- 



Since the jamiton ODE (24) cannot be integrated beyond vm, there is a maximal jamiton 
(for this vs) that connects the specific volumes vm and vr, as depicted by the blue curve in 
Fig. [6j Note that this maximal jamiton is only of theoretical relevance because it is infinitely 
long. However, the above construction gives rise to a one-parameter family (for the specific 
Vs at hand) of jamitons of finite length, whose smooth parts are comprised of the smooth 
part of the maximal jamiton, but which differ in the states that are connected by a shock. 
Specifically, the construction continues as follows: 



7. Select a value r £ (rmin, '''maxy 

8. Determine v"*" < vs and v~ > vs, such that r(v~) = r = r(v+) (see Fig. |5 

9. Calculate desired properties of the jamiton, such as: 



• Solve the jamiton ODE (24) from v+ to v to plot the jamiton in the x variable. 

• Transform the jamiton into the Eulerian self-similar variable rj, by using the relation 

dr] = ^{dx — s dt) ^ ^{dx — udt + mv dt) — ^ (v do" -|- mv dt) = v dx • 
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Figure [g] shows the shape of a jamiton as a function v(ry). 
Compute its length in the rj variable 

L = r]{v-)-r]{v+)= [ i]'{v)dv=[ dv (28) 



w{v) 

and the number of vehicles in the jamiton (in the r/ variable) 

N= r p{v)dr]= r ^dv= r ^dv (29) 

by evaluating a definite integral. Of course, in line with having a complete contin- 
uum description, N is not restricted to integer values. Moreover, in the physical 
road coordinate x, the jamiton length is tL, and the number of vehicles is tN. 

As depicted through the dashed function in Fig. [6j a single jamiton can be continued into 
an infinite periodic chain of jamitons, all of which are connected through shocks. In the r-v 
diagram shown in Fig. [sj this means that one moves continuously along the function r(v) 
from v"*" to v~, and then jumps back to v"^ via a shock. 

In ZND detonation theory, the constant line whose intersection with the graph of r(v) 
determines the two states v"*" and v~, is called Rayleigh line. As one can see in Fig. § the 
function r(v) is convex, and the sonic value vg is the unique point where its slope matches 
the slope of the Rayleigh line. The construction presented here leads us to the following 

Theorem 2. For the PW model, the existence of jamitons with a sonic value vs is equivalent 
(with the exception of the degenerate case, see Rem. [ip to the instability of the base state 
solution v{x,t) = u{x,t) = U{vs)- 

Proof. As argued in Sect. |4.1.1[ a n ecessary condition for the existence of jamitons is that 
U'(vs) > where m = y^— p'(vs). This condition is the complement of the stability con- 



dition (17). Moreover, the construction in Sect. 4.1.2 shows that for every linearly unstable 



value Vs, one can in fact construct jamitons that have this sonic value vg. □ 

4.2. Inhomogeneous Aw-Rascle-Zhang Model. The construction of jamitons for the 
ARZ model is very similar to the construction for the PW model. Again, we seek for solutions 



in the self-similar variable x = . In Lagrangian variables, the ARZ model (13) differs 



from the PW model (18]) only in the fact that {p{v))a is replaced by {h{v))t. Therefore, 



equation (21) carries over, and thus relation (23) holds as well. Equation (22) is replaced by 



r^u'ix) + /i'(v(x))^v'(x) = ^ {U{v{x)) - u{x)) , (30) 

and hence the jamiton ODE ( |24[ ) carries over, with the only modification that the function 
r(v) now reads as 

r(v) = mh{w) + m^v . (31) 



Thus, the condition for the vanishing of the denominator in (24) is now /i'(vs) = — m, and 
therefore for a given sonic value vs, we obtain 

m = — /i'(vs) and s = U{vs) — w-vs • (32) 



For shocks in the ARZ model, the Rankine-Hugoniot conditions (14) yield 

m [v] + [/i(v)] = . 



Thus, with r(v) defined as in (31), the same condition r(v ) = r{v~^) as in the PW model 



applies. The arguments involving the Lax entropy condition transfer exactly from the PW 
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model, as do the considerations regarding roots of 'w{v), as well as the jamiton construction 
principle. We can therefore formulate 

Theorem 3. For the ARZ model, the existence of jamitons with a sonic value vg is equivalent 
(with the exception of the degenerate case, see Rem. [ip to the instability of the base state 
solution v{x,t) = V5, u{x,t) = U{vs)- 

Proof. The condition for the existence of jamitons is that U'(ys) > "t-; where m = — /i'(vs). 



This condition is the complement of the stability condition (19). The remaining arguments 



are the same as in the proof of Thm. [2] □ 

Note that for the ARZ model, the connection between the existence of traveling wave solu- 
tions and the instability of base states was observed by Greenberg [19], albeit without drawing 
connections to ZND theory, using the SCC, or exploiting the geometric interpretations thereof. 



5. Jamitons in the Fundamental Diagram 

We now use the explicit constructions of jamiton solutions of the PW model and the ARZ 
model, presented in Sect. [4], to generate set-valued fundamental diagrams of traffic flow. 
To do so, we consider the solutions of Q or (10), when initializing the model with constant 



base state initial conditions (p, U{p)) plus a small perturbation. For each value of p, for which 
condition ( [l7| ) or ([19]), respectively, is satisfied, the solution evolves towards the constant state 
with density p. In contrast, for each value of p, for which the SCC is violated, the second 
order model (Q or (10), respectively, converges to a jamiton-dominated state. As explained in 



Sect.[3j the former statement follows from known rigorous results for second order traffic flow 
models with relaxation terms, while the latter statement is based on numerical experiments. 

In order to construct a fundamental diagram, we start with a plot of the equilibrium curve 
of traffic flow 

Qeqip) = pU{p) 

in the flow rate vs. density plane. For every value of p that satisfies the SCC, we can record 
a point {p,Qeq{p)) on this equilibrium curve, representing the stable base state solution. In 
turn, for every value of p that violates the SCC, we sketch a jamiton line segment, which is 
defined as follows. 

As derived in Sect. |4j a jamiton is non-constant in p and u, and therefore it does not 



represent a single point in the flow rate vs. density plane. Instead, due to relation (23), which 
implies that u = ^ + s, each jamiton is given by a segment of the line 

Q = m + sp 



in the flow rate vs. density plane. An example of such a jamiton line is depicted in Fig. 19 , Due 



to the Chapman- Jouguet condition (26), every jamiton goes through the equilibrium curve 
Qeqip) at the sonic density ps- For each choice of sonic density ps (that violates the SCC), 
the corresponding maximal jamiton is represented by a line segment between (pM, + •sPm) 
and {pB.,m + spR,), where m, s, pu = I/vm, and = 1/vr are constructed in the ffist 6 



steps presented in Sect. 4.1.2 Since w{vm_) = 0, where w{v) is defined by ( p5| ), it follows that 
PmU{pm) = rn + spu, i-e., the point {pM,m + spu) lies on the equilibrium curve Qeqip) as 
well. Because Qeq{p) is assumed concave, the maximal jamiton line therefore lies below the 
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equilibrium curve for densities pM < p < PS; ^'iid above the equilibrium curve for densities 
PS < P< PR- 

Moreover, the slope of each jamiton line, s, equals the propagation speed of the jamiton in 
a stationary frame of reference. In particular, if s > 0, the jamiton is moving forward, and if 
s < 0, the jamiton is moving backward, relative to a stationary observer. Here, the former 
case is associated with light traffic and the latter case with heavy traffic, as formalized by the 
following theorem, which provides information about the structural properties of a jamiton 
FD. 

Theorem 4. For both the PW and the ARZ model, the following two properties hold: 
(i) The jamiton velocity s decreases with the sonic density ps- 

(a) At the transition points between stability and jamiton solutions, the jamiton lines are 
parallel to the equilibrium curve. 



Proof. For (i), we show equivalently that s increases with the sonic specific volume vs. Due 
to Thm. [2] and Thm. [sj jamitons exist when U'{vs) — m{vs) > 0. Due to (27) and (32), we 
have that s(vs) = U{vs) — '"T'(vs)vs. Hence 

s'(vs) = f^'(vs) - m{vs) - m'(vs)vs > -m'(vs)vs . 

Moreover, since we assume that p"{v) > for the PW model, and h"(v) > for the ARZ 
model, we find that m'(vs) < and therefore s'(vs) > 0. 

To show (ii), we note that due to Thm. [2] and Thm.jsj the boundary between stability and 
jamitons occurs, in Eulerian variables, when —U'{ps)p^= m(ps). Hence 

s{ps) = Uips) - ^ = Uips) + U'{ps)ps = Q'iPs) , 
which proves the claim. □ 



5.1. List of Examples. Before presenting the construction of jamiton FDs (and extensions 
thereof, see below) in detail, we first collect the parameters of the examples that we shall use 
as test cases. In all examples, we consider pmax = l/7.5m and Umax = 20m/s, except when 
explicitly noted otherwise. Specifically, we consider the following four scenarios: 

PWl The PW model with a quadratic equilibrium curve, i.e., a linear desired velocity 
function ([2]), and a singular pressure of the form (18), with (3 = 4.8m/s^. This model 
is used in Figs. [7 10 



PW2 The PW model with a non-quadratic concave equilibrium curve of the form 

Q{P) = C (5(0) + (5(1) - 5(0))p/Pmax - gip/Pmax)) , 



(33) 



where the choice of the function g{y) = Wl + (^^)^ has the effect that (33) can be 



interpreted as a smooth version of the piecewise-linear Newell-Daganzo flux [^[ [8]. 
The parameters_ are chosen such that a reasonable fit with the measurement data, 

the 



shown in Fig. [ij is achieved: c = 0.078/9 



max^max; 



and A 



logarithmic pressure (18) is chosen, this time with /3 
Figs. 



8m/s2 



Jo- Again, 
This model is used in 



11 14 



ARZl The ARZ model with the same equilibrium curve (33), and with the singular hesitation 
function (20), which was suggested in |3], with 7=5 and /? = 8m/s^. This model is 



used in Figs. 15 - 18 
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ARZ2 The ARZ model with the equihbrium curve (33), and with the singular hesitation 
function 



h{p) = /3 



\7i 



(1 - P/Pn 



^72 



where 71 = 72 = and /3 = 8m/s^. This generalization of (20) yields a better 
qualitative fit of the unstable regime with measurement data (see Sect. 5.5). This 



model is used in Fig. 21 



Note that some of the model parameters (pmax) ^max) StiQ chosen to be within ranges of values 
that are realistic for real traffic, some (c, 6, A) to match specific data, and some (/3, 7^) to 
yield reasonably sized jamiton regions. 



5.2. Maximal Jamiton Envelopes. The relations described in the beginning of this section 
give rise to the following construction of a set-valued fundamental diagram for the PW or the 
ARZ model, defined by a desired velocity function U {p) and the traffic pressure p{p) or the 
hesitation function respectively. 

Loop over all densities < p < /Omax- For each density p: 

1. Check whether the SCC is satisfied. If it is, draw the equilibrium point {p,Qeq{p)). 



2. If it is not, compute m, s, pu, and /Or (see Sect. 4.1.2), and plot a straight line between 
the points {pu, m + spy\) and (pR, m + spji). 

An example of such a jamiton fundamental diagram is shown in Fig. [Tj It consists of segments 
of the equilibrium curve (wherever p satisfies the SCC; shown in black), and a set-valued region 
in the flow rate vs. density plane, in which the equilibrium curve lies (for all p that violate the 
SCC). We call this set- valued region the jamiton region, which is formed by a continuum of 
jamiton line segments, shown in blue in all figures that present a jamiton FD. The boundaries 
of the jamiton region are formed by an upper envelope curve above the equilibrium curve and 
a lower envelope curve below the equilibrium curve, both of which are depicted in magenta. 

Theorem |4| together with the concavity of Q{p), implies that different jamiton lines cannot 
intersect above the equilibrium curve. Hence, the upper envelope of the jamiton region is a 
curve that is formed by the points (/3r, m + sp^), where m, s, and />r are parameterized by 
ps- In contrast, different jamiton lines do intersect below the equilibrium curve. Therefore, 
the lower envelope of the jamiton region is formed by intersection points of jamiton lines. In 
Eulerian variables, this curve, wherever it is below the equilibrium curve, is given by {p* , Q*), 
where 

P*{ps) = and Q*{ps) = m(/Js) + s(ps)p*(ps) 

s'ips) 

represent the parametrization by ps- 

This construction of a maximal jamiton FD is applied to the test cases PWl, PW2, and 
ARl, described in Sect. 5.1, The results are shown in Figs. [7j [TTj and 15 One can clearly 



observe the structural properties described above. In comparison with the FD obtained from 
measurement data (see Fig. [T]), one particular disagreement is apparent: in the jamiton FD, 
the upper envelope curve extends much further than in the data FD. This disagreement can 
both be explained and remedied, as done below. 
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Figure 7. Maximal jamiton FD (i.e., no tem- 
poral aggregation) for ttie example PWl (see 
Sect. 



0.4 0.6 
density p / p 

Figure 8. Jamiton FD with short aggregation 
time At — T for the example PWl (see Sect. 5.1 ). 
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Figure 9. Jamiton FD with long aggregation 
time At — 8t for the example PWl (see 
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Figure 10. Jamiton FD with effective flow rate 
(i.e., aggregation over complete jamitons) for the 



Sect. 5.1). 



example PWl (see Sect. 5.1). 



5.3. Emulating Sensor Measurements. For each jamiton line in the fundamental dia- 
gram, the segment below the equilibrium curve corresponds to the (possibly infinitely) long 
low density piece of a jamiton, while the segment above the equilibrium curve corresponds to 
a comparatively short piece of the jamiton, at which the density is high (see Fig. [6]). In addi- 
tion, as one can see in the example jamiton shown in Fig. [6| the high density part frequently 
comes also with a large negative density gradient along the jamiton curve. 

This can explain the apparent qualitative disagreement between the upper envelope curves 
in the jamiton FD (see Figs. [7 11, and 15) and the data FD (see Fig. [T]). Measurement data 
used to generate fundamental diagrams are typically obtained by stationary sensors on the 
road that measure the average density and average flow rate during a time interval At. In 
other words, each point in a data FD corresponds to an aggregation of the traffic situation 
over a non-zero time interval. In order to obtain reliable estimates for the flow rate (which 
is obtained by counting vehicles), the aggregation time At must not be too small. A typical 
value is At = 30s, as is used for instance in the Mn/DOT RTMC data set [H], as described in 
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Figure 11. Maximal jamiton FD (i.e., no tem- 
poral aggregation) for the example PW2 (see 
Sect. 
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Figure 12. Jamiton FD with short aggregation 
time At — T for the example PW2 (see Sect. 5.1 ). 
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Figure 13. Jamiton FD with long aggregation 
time At = lOr for the example PW2 (see 
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Figure 14. Jamiton FD with effective flow rate 
(i.e., aggregation over complete jamitons) for the 



Sect. 5.1). 



example PW2 (see Sect. 5.1). 



Sect. 1.2 Hence, in order to obtain a better agreement with sensor measurements, we must 
represent the jamiton FD in a way that is more faithful to the actual measurement process. 
We therefore perform an averaging over jamiton solutions, as follows. 

For a given sonic density ps, we consider a one-parameter family of traveling wave solutions, 
each being an infinite chain of identical jamitons, where the parameter is the length of each 
individual jamiton, measured from shock to shock. For each such periodic traveling wave 
solution, we consider temporal averages of duration At at a fixed position x on the road. 
Because the chain of jamitons moves with velocity s, these temporal averages can be translated 
into averages over the Eulerian similarity variable t] = Specifically, the average density 
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Figure 15. Maximal jamiton FD (i.e., no tem- 
poral aggregation) for the example ARZl (see 
Sect. 
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Figure 17. Jamiton FD with long aggregation 
time At = lOr for the example ARZl (see 



0.4 0.6 
density p / p 

Figure 16. Jamiton FD with short aggrega- 
tion time At = T for the example ARZl (see 
Sect. 
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Figure 18. Jamiton FD with effective flow rate 
(i.e., aggregation over complete jamitons) for the 



Sect. 5.1). 



example ARZl (see Sect. 5.1) 



IS 



1 / n( 7^ 4\ A+ 1 / „/' X — St^ 



(34) 



rv 

Jri-Ari 



where Arj = is the averaging "length" (in the r] variable), and rj = - parameterizes the 
"position" (in the rj variable) of the averaging interval. Similarly, the average flow rate equals 



rv rv 

Q = ^ Qiv) ^V = ^ m + spiv) dr] = m + sp , 

Jr^—Ari Jri—Arj 



(35) 
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which imphes that any average pair {p, Q) lies on the corresponding jamiton hne Q = m + sp. 
Note that these averages depend only on the dimensionless quantity 

At 



a = — 
r 



rather than At or r alone, which measures the temporal averaging duration in units of the 
drivers' relaxation time r. 



For sufficiently large At and for sufficiently short jamitons, the averages (34) and (35) 
may involve one or more shocks. In this case, it is helpful to split the integrals into multiple 
segments, each of which corresponds to a smooth jamiton curve. For each of these segments. 



it follows from (29) that 



rj pv{ri) 



r'(v) 



P=^ Pi^l) dr? = ^ / dv , 

which can be computed efficiently using composite Gaussian quadrature rules. 

Using these averaging procedures, we can now, for each choice of a = — , produce an 
aggregated jamiton fundamental diagram, by plotting on each jamiton line segment all possible 
pairs {p,Q) that can arise as averages. As a consequence, the jamiton line segment above 
the equilibrium curve is shortened, since the peak density /jr is averaged with lower density 
values. In turn, the line segment below the equilibrium curve is not shortened: due to the 
infinite length of the maximal jamiton towards pM, averages can always become as low as pM 
itself. 

Aggregated jamiton FDs are shown for the example PWl in Figs. [8] and [9j for the example 



PW2 in Figs. [12] and 13, and for the example ARl in Figs. 16 and|17[ In all three examples, 
the ffi'st figure corresponds to a short aggregation duration of At = r, and the second figure 
corresponds to a long aggregation duration of At = 8r or At = lOr, respectively. The 



maximal jamiton FD, constructed in Sect. 5.2, is recovered in the limit a — )• 0, for which the 



averaging length goes to zero, and thus averages turn into point evaluations, i.e. 



V 



A comparison of the referenced aggregated jamiton FDs with the data FD (shown in Fig. [T]) 
shows a better qualitative agreement than in the case without the aggregation. In particular, 
the very high values of Q, corresponding to high-density parts of jamitons, are significantly 
reduced. An interesting feature that can be seen in Figs. [9, 13, and 17 is the cusp at the 



PR point of the jamiton line that is parallel to the p axis. The reason for this is that for 
s = 0, the jamiton is not moving along the road, and thus aggregated values are identical to 
un-averaged point evaluations. 



5.4. Effective Flow Rate. In all jamiton FDs shown in Figs.[^-[9| Figs. [TT] - [T3l and Figs. [15]- 



17, the jamiton region lies partially above the equilibrium curve. This is due to the fact that 
on their upstream side, jamitons possess densities that exceed the sonic density (see Fig. [6]). 
Therefore, if averages involve more of this high-density part than of low-density parts, then 
the aggregated density and flow rate pair (p, Q) may lie above the equilibrium curve. 

To obtain the true effective flow rate that a chain of jamitons yields, one needs to perform 
the aggregation over inflnitely long time intervals (we are neglecting the singularity at s = 
here), or equivalently, over full jamitons, i.e., from shock to shock. This allows one to 
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Figure 19. The maximal jamiton (blue line) and Figure 20. Secants (orange lines) to the func- 

an actual jamiton (brown line) in the Q-p plane. tion w{\) in the u-v plane, connecting each of 

Both jamitons lie on the same straight line that the jamiton volumes v+ and with the sonic 

goes through the sonic point (ps, Q{ps))- volume vg. 



investigate the important questions: 1) Do jamitons always lead to a reduction in the effective 
flow rate, compared to uniform traffic flow with the same average density? 2) Can one 
construct physically realistic examples in which jamitons actually increase the effective flow 
rate? 

To answer these questions, we flrst apply the full-jamiton averaging to the three test cases 



presented in Sect. 5.1 The resulting effective flow rate jamiton fundamental diagrams are 



shown in Figs. [Toj [T4j and 18 One can observe that in all three examples the jamiton region 



below the equilibrium curve remains unchanged, while it has fully vanished from above the 
equilibrium curve. The former observation is a consequence of the fact that on the one hand, 
very short jamitons possess average densities very close to their sonic density, and that on the 
other hand, very long jamitons possess average densities close to their minimum density. In 
contrast, the latter conclusion, which means that in the considered examples jamitons always 
reduce the effective flow rate, does not follow immediately from the shape of the jamitons. 
However, the theoretical answer to this question is provided by the following 

Theorem 5. The effective flow rate of a periodic chain of identical jamitons is always lower 
than the flow rate of a uniform base state solution with the same average density. 

Proof. The average flow rate Q being less than the flow rate of the uniform flow with the 
same average density, Q{p), means that the averaged pair {p,Q) lies below the equilibrium 
curve (see Fig. [19]). Thus, we need to show that averages over full jamitons satisfy p < ps, 



where p = N/L, and and L are given by (28) and (29). 



Due to the Rankine-Hugoniot shock conditions, given in Sect. 4.1.1 , a jamiton lies between 



two states that possess the same value of r(v) (see Fig. [s]). The function r(v), given by (25) 
or (31), respectively, can be inverted on two branches: let v~^{r) denote its inverse function 
on [vm,vs], and v~(r) its inverse on [vs,vp{]. We can therefore parameterize the family of 
jamitons by f G [^'(vs), "''(vr)]. For each jamiton, we obtain 

Pif) = N{r)/L{r) , where N{r) = / ^ dv and L{f) = / v^ dv , 

Jv+{r) Jv+{r) 
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where the jamiton hes between the specific volumes v~^{r) and v (r). 

We need to show that N{r)/L{r) < ps, or L{r)/N{r) > vs, which is equivalent to the 
condition 



(v- vs)^^dv > . 

v+(r) W{v) 

After changing the variable of integration to f = r(v), this becomes 

v~(r)-vs v+(f)-vs 



/ / df > , (36) 

r(vs) w{v (r)) w;(v+(r)) 
where the two quotients arise from the two branches of the inverse of r(v). We now show that 



the integrand in (36) is positive, which proves the claim. 



As depicted in Fig. 20, consider two secants of the function w{y): the first goes through the 



points (v+, u)(v"'")) and (vs, 0), and the second goes through the points (vs, 0) and (v , w{v )) 



(omitting the argument f for notational economy). By assumption (see Sect. 2.3), U'{v) > 0, 
thus both secant slopes are positive. Moreover, U"{v) < 0, thus w{v) is concave, and therefore, 
the slope of the first secant is larger than the slope of the second, i.e. 

w{v~) ^ tt;(v+) 



V — Vs V""" — Vs 

Observing that the reciprocals of these slopes are the same as the two terms in the integrand 



of (36) concludes the argument. □ 



5.5. Mimicking Real Sensor Measurement Data. To conclude the series of examples, 
we now consider the ARZ model with parameters selected such that its aggregated jamiton 
FD resembles the measurement data in Fig. [l] and the phase diagram in Fig. [2| at least 



qualitatively. We do so by using example ARZ2, as described in Sect. 5.1, with the following 
two modifications: Umax = 19.2m/s, /3 = 12m/s^. In addition, we choose a = 6, i.e., for a 
data aggregation time of At = 30s, this corresponds to a drivers' relaxation time of r = 5s. 



The results of these choices can be seen in Fig. 21, which shows the aggregated jamiton FD 



together with the measurement data. A strong qualitative agreement is apparent. 
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6. Conclusions and Outlook 

We have demonstrated that both the classical Payne- Whitham model and the inhomoge- 
neous Aw-Rascle-Zhang model of traffic flow possess an intrinsic phase transition behavior 
between stable and jamiton-dominated regimes. This structure gives rise to set-valued fun- 
damental diagrams, and we have provided an almost fully explicit way to construct such set- 
valued fundamental diagrams from the parameters governing the second order traffic model. 
The set- valued part of the fundamental diagram represents families of traveling wave solutions, 
so-called jamitons. Various ways to perform averaging on these solutions (thus emulating sen- 
sor measurements) lead to different types of envelopes for the set-valued region. Appropriate 
choices for the model parameter functions and for the type of averaging lead to fundamental 
diagrams that, at least qualitatively, agree with fundamental diagrams obtained from sensor 
measurements. 

The phase transitions between stable and jamiton-dominated traffic flow are related to the 
sub-characteristic condition, that was introduced by Whitham in 1959, and the construction 
of jamiton solutions is based on the Zel'dovich-von Neumann-Doring (ZND) theory of reacting 
gases. Through both relations, interesting connections between traffic models and many other 
physical phenomena can be drawn (see [15] for examples). 

Having a systematic approach for the construction of the envelopes of a set-valued funda- 
mental diagram from a given traffic model could also help in the inverse problem of deter- 
mining desired velocity functions and traffic pressures from observations and data, as follows. 
While (at least for congested traffic) these parameter functions cannot easily be measured 
directly, envelopes of data points (modulo outliers) can certainly be obtained directly from 
measurements. Thus, the relationship between the fundamental diagram and jamiton solu- 
tions established in this paper could serve as means to re-construct parameter functions so 
that the model reproduces the observed envelopes. 

The methodology presented here allows the complete construction of fundamental diagrams 
from second order traffic models. Yet, certain aspects require a deeper investigation. Most 
prominently, the dynamic stability of the jamiton solutions is not considered. While numerical 
studies (e.g., [IHIIIS]) suggest that many jamiton solutions are in fact stable attractors, this 
observation may not generalize to all jamiton solutions. The constructions conducted in 



Sect. 4.1.2 can, in principle, generate jamitons of arbitrarily short, and of arbitrarily long, 
wavelengths. Both are most likely dynamically unstable: a chain of very short jamitons 
is essentially a uniform density state plus some small sawtooth perturbation; and a very 
long jamiton possesses a very long "tail" of almost constant density. In both cases, if the 
corresponding density is itself unstable, small perturbations will amplify and thus destroy the 
structure. While it remains a topic of future research to analyze the structural stability of 
jamiton solutions, it can certainly be concluded that the removal of long wavelengths might 
result in interesting modifications of the resulting fundamental diagrams. 

Systematic numerical simulations can be employed to study the dynamic stability of jami- 
ton solutions. In this fashion, fundamental diagrams can be systematically generated through 
numerical simulations (as conducted for instance in [47j ) . Therefore, numerical studies could 
provide a missing link between fundamental diagrams obtained from sensor measurements 
and those constructed via jamiton solutions. 

One limitation of this study is the restriction to the Payne- Whitham model and the Aw- 
Rascle-Zhang model only. Many more general types of second order traffic models do not fall 
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into either of these two categories, in particular models that are partially based on measure- 
ment data (e.g., generalized Aw-Rascle-Zhang models [13j). In a forthcoming paper [16j, we 
plan to extend these results to a more general class of second order traffic models that contain 
both the PW pressure and the ARZ hesitation function at the same time. 
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